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Abstract. 

The question of allowed signs of the static dilectric function for exactly 
homogeneous ground states in many body systems is further analyzed. The 
discussion is restricted to zero temperature situation. Firstly, it is argued 
that the positivity is equivalent to the requirement that the work needed to 
adiabatically connect an static charge density, be positive. In particular the 
rule seems to be a fully appropiate stability condition if the system interact 
with a freely moving but Coulomb interacting compensating jellium. This 
situation is common in thermodynamical equiliubrium discussions, in which 
the system can be stable or not in dependence of the relaxation of some 
external constraints. Further, an argument based in a quantized electromag- 
netic field treatment is also given. It requires e > as an intrinsic stabil- 
ity condition for the charged system independently of the jellium dynamics. 
Possible implications for the more complex finite temperature situations are 
commented. 
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1 Introduction. 



The problem of the positivity of the static dielectric function have been 
debated in the literature in close connection with a related but not equivalent 
question about the validity of the dispersion relations [1-4]. In the works [3] 
and [4] it has been argued that negative values of e(k) are allowed by stable 
translation invariant physical systems. The conclusion was extracted from a 
discussion of appropriate effective actions which should have minimum values 
at equilibrium. 

In a previous work [5] we have considered the requirement for a posi- 
tive work in adiabatically creating a fluctuation of the charge density in the 
background. The conclusion was that such a condition implies the positivity 
of e(k). Then, in this work we intend to clarify the connections among the 
results of [3], [4] and our work [5]. Only the null temperature situation will be 
considered here. The harder physical problematic at T ^ will be discussed 
elsewhere. 

The work start in Section 1, from adapting the variational derivation of 
the positivity of the effective action, to the problem [6]. The considered 
physical system consists of several Coulomb interacting species of charged 
particles in interaction with a compensating jellium. Then, the adiabatic 
connection of the auxiliary electrostatic external source allow to obtain the 
same stability conditions of [3] and [4] from the positivity of the effective 
action. The variational discussion being used clarifies the meaning of the 
positivity for the work done in connecting the charge functuation [5]. It is 
simple equivalent to requiring that the energy E[p] of the state after the 
adiabatic connection of the sources becomes greater than the ground state 
energy E g . The inequality e(k) > directly follows from this requirement. 

A first conclusion then follows. It refers to systems in which the jellium 
charge density is freely moving but Coulomb interacting with the particles, 
that is, systems in which the jellium homogeneity is assured self consistently 
by the Coulomb interaction. In such situation if E[p] can be lower than the 
ground state energy E g , the system can relax by itself to an inhomogeneous 
configuration by generating a jellium density modulation. Therefore, in order 
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that the whole system: particles + jellium, be stable the positivity of e(k) is 
needed. 

Next, in Section 2, after considering the electromagnetic field as quan- 
tized, the positivity of e(k) follows as an intrinsic stability requirement com- 
ing from the internal dynamics of the charged particle system and being 
independent of the jellium nature. 

In the summary, the main conclusions are reviewed and comments on the 
posible implications for the finite temperature case are given. 



— * 

2 e(k) > and positive work conditions 

Let us consider a many body system consisting of r Coulomb interacting 
species of charges. Its quantization is described by the transition amplitude 

Z[<f>] = 7|| =< n a> H 5 ^' out ><P (!) 

where \Q g > is the ground state associated to the Hamiltonian in the absence 
of the external electrostatic field and it is assumed that the static field (f>(x) 
is adiabatically connected from \Q g ,in > at the begining, further maintained 
during a long interval T and after disconnected when the state becomes 
\Q g , out > [6]. 

The hamiltonian after introducing the field <f>(x, t) is given by 



J 8=1 J 8=1 



+ 



In 
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:,t=l 



e\x — y\ 



(2) 



where means normal ordering, Q s are the kinetic one particle operators 
for each componet and po a possibly needed jellium background density. The 
external field 4>(x) is assumed being generated by an electrostatic charge 
density through the Poisson equation 
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/4-ir 
T^^p(y)d\ (3) 
\x y\ 

Let us consider a translationally invariant ground state of the system 
\Q g > and the adiabatic connection of a specific density p(x) under the above 
defined scheme. Then, the ground state to ground state transition amplitude 
given by Z in (1) should obey 

< Q g ,in\n g ,out >= exp ( - ^-A£[0]t), (4) 

with AE[(f)] = E[4>] — E g and E[<p] and E g are the ground state energies 
in presence and absence of the source p respectively. For the generator of 
connected Green functions W it follows at large L values 

W[<f>(x)] = -ih\nZ[(f)(x)} = -(E[(f)} - E g )T. (5) 

Next, consider the search for the minimum of the mean value of the energy 
in the state space under the restriction that the mean value of the internally 
produced charge density takes a fixed static value pi(x) that is 

<fl|Eg Qs^ti^sjx) — po I O > _ 

The usual normalization condition < Q\Q >= 1 will be also required. 

Alternatively the extremum to be found can be obtained by using the 
method of Lagrange multipliers and minimizing the quantity 



< n\H\n > -a < n\n > - / d 3 xp(x) < n\J2qsi>f(x)i>s(x)-po\n >, (7) 

J 8=1 

in which a and (3{x) are to be determined from imposing the normalization 
and mean field condition (6). The extremal of these expressions over arbitrary 
states leads to the following Schodinger equation for the solution \Q >. 

(H - [ d 3 x(3(x) j2Mt(x)M%) - P.)) P >= a\n > . (8) 
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But, on another hand the state which is obtained after the adiabatic 
connection of a density p(x) is also abeying the Schrodinger equation 

(H-J d 3 x^x)(Y,^t(^s(x) - po)|fy >= E[<t>]\% > (9) 

s=l 

Therefore, if the auxiliary charge p(x) is taken sufficiently small in order 
that no crossing of energy levels occurs in connecting it, the state > is also 
the ground state in presence of the source p(x). Thus, minimal conditions 
under the constraints (6) and the normalization condition are obeyed by the 
state \ flfj) > if the multiplier values are selected to be 



a = E[<P Pi ] 
(3{x) = (10) 

where <p Pi (x) is such a field that forces the internal density 

s=l 

to be equal to Pi(x). Here the central point comes. Taking the scalar product 
of (9) with > the minimal values for the mean energy, once the field 

4>{x) is fixed, takes the form 



= AE[(/>] + J d 3 xp t (x)(j)(x) 

= ~{W[<f>\- j d 3 xp t (x)m) + E g 

= E g -^V[pl (11) 

where <f>(x) = <f> Pi {x) and the effective action T[pi] as a functional of any space 
time dependent field Pi(x), x = (x, t) is given by the Legendre transformation 
from the variables <j)(x) to the mean fields Pi(x) 

T[p t ] = W[<t>{x,t)\ - J dx Pi {x,t)4>{x,t) (12) 
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pl[x) ~ S(j>{ x y 

, 5T\Pi] | ^[0] | __,(4)/ _ x 



But, < iJ >q in (11) is the minimum value of the mean energy under a 
restriction, therefore it should be necessarely greater than the ground state 
energy which is the fully unrestricted extremal value. Thus, the stability 
condition for the assumed translational invariant ground state becomes 

< H > n0 -E g = ~T[ Pi (x)\ > 0, (13) 

which requires the positivity of the effective potential under variation of the 
internal density Pi(x) 



V[pi(x)] = - lim 

1 — >oo 



(14) 



Thus, the functional formed by the second derivatives of V should be a 
positive definite kernel. This is the stability condition obtained by Kirshnitz 
[3] and it should be certainly obeyed. As it was discussed in [4] condition 
(14) also allows for negative values of the dielectric function. It will be seen 
in what follows. Integrating the last of relation (12) at vanishing = Pi = 
over both time t x and t y , dividing by the integration interval T and making 
use of the assumed translation invariance, it follows for the second derivative 
of V[ Pi \ 

jd 3 z(j dtY {2 \x - z,t)) J dt'W {2 \z - y,t') = -5 (3) (x-y), 

Jt'sSm^™™*-''"- ^*-*' (15) 



in which 



^ 5 2 W 
r (2) (x-y)= I U =0 , (16) 

dp i {x)5p i (y) 
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= (17) 

and W^ 2 \x, u) is the temporal Fourier transform of W^ 2 \x). 

y That is, the second derivative kernel of the effective potential is the in- 
verse of the temporal Fourier transformation of W^ 2 \ The Fourier transform 
can be evaluated as usual by introducing the completion formula of the state 
space between the two operators in the relevant time ordering. It produces 



WW(x-gO) = ^2Re[<0\j2Qs^f(x)M^)\n> 

1 



^ n^O s=l 



( e n _ e o) 



<n\J2Qs^f(y)My)\o>], (is) 



s=l 



where |0 >= > and e = E g . 

For the space Fourier transform of (18) ot follows 



»*'M)4E* (19) 

with Pi(k) equal to the total density operator Y, s Qs^fffl^six)- This expres- 
sion for shows that the kernel W^ 2 \x — y) is positive definite whenever 
Eg is really the ground state energy, that is E g < e n . Then, the kernel formed 
the by the second derivatives of the effective potential should be also positive 
definite 
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d*xd 3 y Pl (x) f^f J Pl (y) > 0, (20) 

for an arbitrary density Pi(x). 

Let us consider now further independent restrictions required by the sta- 
bility of certain systems. For this purpose it will be made explicit below 
how the source Pi(x) can be considered as a fluctuation of the jellium charge 
density po- Note that the fixing of a definite jellium density is an external 
constraint on the system. In most of the real systems the charge density of 
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the homogeneous jellium is fixed self consistenty by the interaction with the 
charged particles. 

Consider the interaction of the system with a modified jellium in the form 

lr r 
I = - d 3 xd 3 y(CZl q s ipf(x)ip s (x) - p - p(x)). 

1 J s=l 

r 

^(^-$(H^ s + (y)^(y) - Po-p{y)), (21) 

s=l 

where the ordering of the density operators is not relevant if the Coulomb 
interaction is regularized through a parameter e in such a way that it vanish 
at x — y = 0. Also it can be assumed that v e remain being positive function 
v e (x — y) > 0. 

Relation (21) can be also written in the form 



f d 3 xd 3 y q s ipf(x)^ s (x)v e (x -y)J2 Qs'&ttififaiy) 



s=l s=l 



- J2^t(^s(x)v e (x-y)(p + p(y)) 

s=l 

+ ^(po + p)(S)v e (x-y)(p + p(x)). (22) 

But due to the constancy of po it follows that on any state with fixed 
number of particles and in a thermodynamical limit sense 

f r 1 f - 1 

/ d 3 xd 3 y Vftf s + (f)i(i)- -p \Q >= / d 3 x Q— p |O > 

J \x-y\ J Fl 

= f d 3 xd 3 yp —^—p \n >, (23) 
J \x — y\ 

where Q is the total charge operator. 

By also considering that the total jellium charge does not vary in the 
fluctuation, that is 

J d 3 xp(x) = 



, it follows 



/ = / d*xd 3 y(± qs ^(mZ) Ve{S 9 ®- T,q.1>tm.W (24) 

J s=l Z s=l 

r 111 

- (^2(ls^f(x)^ s (x) - p )(f)(x) + —d(f)(x)d<j)(x) + -p — — -p ). 

s =i 47r * \ x — y\ 

Then, the Coulomb interaction operator of the system in (2) coincides 
with I apart from the necessary divergent self energy of the jellium plus the 
electrostatic self energy of the fluctuation. Clearly, the absence of the last 
term in (24) means that it is added in the energy E[<j>] in (9). 

Finally let us consider the work done by the external systems in order to 
impose a constraint consisting in the variation p in the jellium density (or its 
equivalent: a fixation of evalue pi(x) of the internal density). 

Note that the work E that should be performed to create, let say the 
mean value of the internal density Pi(x) should be equal to the energy of the 
state after connecting (p Pi (x) plus the energy of the electrostatic field (f> Pi (x) 
minus the ground state energy E g . That is 

S[0 p J = E[<t> Pi ] + 1-J d'xd^d^Xx) - E g . (25) 

Note, that this quantity differes radically from the effective potential (14). 
From (5) it follows 

5E(4>) 
5(f) Pi (x)5(f) Pi (y) 
W {2) (x-y,tu = 0) 

Re[< 0| E s q s ^ + {x)^{x )\n >< n\J2t Qtipt ~{y)My)\® >] 

(26) 

Hence, the second derivatives of £ are given by 



£ (2) [0p 



\ " o 



Tin 
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= -±V^(f- S) -ig2*[<0|i:rf(fW,(f)|n> 

— < n\Y,q t i,tmM\0 >]• (27) 
Further, its spatial Fourier transform takes the form 



eA; 2 



4vr e (jfe) : 



(28) 



in which as, usual for the dielectric e, reponse R and polarization P functions 
the following relations take place [7] 



P(k) 



l + v(k)R(k) = 



e(jfe) ~ ' l-v(k)P(k)' 
R(k) 



1 + v(k)R(k) ' 



i2(jfe) = P ffl - (29) 

1 - v{k)P{k) 



2 y-, | < 0\f>k\n > 



|2 



Notice that the stability condition following from the positivity of the 
effective action is given by [3] 

-i?$= * -L( C (£)-1)>0 (30) 
v(k) e(k) 

which implies the inequalities 

e(jfe) > 1 or e(jfe) < 0. (31) 
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Then, we arrived to the main point; the requirement that a positive 
work is always needed to produce a definite value of the internal charge 
density Pi(k) leads to the positivity of of the dielectric function e(k) or its 
inverse. However, this condition should not be directly interpletable as a 
stability one. In order to be so, the jellium should be a freely moving system 
which homogeneity is self consistently determined by its strong Coulomb 
interaction with the particle components. Although, such should be the case 
in many systems, it is not the general one. For example, if the charges of 
the components completely cancels, then p might describe only an external 
auxiliary charge. Then, the work to stablish it is not excluded that could 
be negative. However, if the jellium is a freely moving and non-vanishing, it 
becames clear that the full system,that is particles + jellium, can diminish 
dinamically its energy by adopting an inhomogeneous configuration whenever 
e(k) can be negative. 

The existence of further intrinsic physical conditions requiring the posi- 

— * 

tivity rule e(k) > will be considered in the next section where the electro- 
magnetic field is considered as a dynamical one. 

3 eik) and dynamical photon fields 

In this section, an additional support for the necessary positivity of the static 
dielectric function is given. As it can be reminded, in last section such a 
condition only followed after assuming that the jellium compensating charges 
were dynamical ones. Then, the stability of the whole system required e(k) > 
0. However, as it will be discussed below , it seems that the internal stability 
of the charged system also requires that rule, no matter the character of 
the jellium dynamics. The argument follows from a treatment in which the 
electromagnetic field is considered as quantized. Also the system will be 
described at zero temperature by the euclidean quantum field theory limit of 
the Matsubara finite temperature theory [8] -[9]. Let us consider the partition 
function of a charged system in interaction with the electromagnetic field as 
given by 



Z[j] = Tr 



exp(-/3H)T (exp J ~ J dx j^x)A^x) 
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— J V^+V^VA^exp S-^-Jd 4 x ]^{x)A^x) 



(33) 



where the euclidean action for the charged system can be considered in the 
form [9] 

s = Y c \l do ° 4 J d3x ( c f^s> ^» A ^ - F ^ 

in which F^ v = d^A v — d u A^, ip s are the fermion fields describing the 
s = l,....r species of charged particles and A^ is the electromagnetic field. 
The operator H is the system hamiltonian in the extended state space in- 
cluding the non-physical photons as well as the auxiliary ghost particles. 
The coordinate X4 is a real one and takes values in the interval (0,/3) with 
f3 = hc/kT. The a-dependent term in (34) is associated to the Fadeev- Popov 
quantization of the gauge electromagnetic field. In the case of the electro- 
magnetic field, the contribution to the partition function (33) of the auxiliary 
ghost particles compensating the scalar and longitudinal unphysical degrees 
of freedom, becomes a field independent determinant. Here it is assumed 
to be included in the normalization constant N in (33). For concreteness, 
the fermion part of the lagrangian can be taken to be the summ of r terms, 
each one corresponding to the interaction of non-relativistic particles with 
the quantized electromagnetic field, that is [10] 



-(X>V', + W</'„W ~Po)A t (x). (35) 

S = l 

Expression (35) for Cf should be very accurate for the description of 
typical condensed matter physical systems in which the dynamics of the 
charged particles can be considered as non-relativistic. 

The mean electromagnetic fields can be obtained through the logarithmic 
derivative of Z with respect to j as 



(AM) = -no stn m - 5 wl3] 



5 jp(x) 5 j^x) 

fVjj+V^VA^ A^x) exp[ff -Jdx j^AJ 
J Vij+VijVA, exp [S-Jdx jpAJ 



(36) 



in which the generating functional of the connected Green functions W\j\ is 
defined by 

W[j] = -he in Z[j] (37) 
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Up to now, the finite temperature situation is yet being considered. Fur- 
ther, the T = limit will be assumed. 

It can be remarked that the variety of subtle points associated with the 
thermodynamical T^O situation need for a more involved discussion in order 
to properly interpret the numerical or experimental data and the existing 
viewpoints [4]- [5]. These question will be considered elsewhere. 

Now, it is possible to perform a Legendre transformation to consider 
the mean fields as the new independent variables in place of the sources j, 
through 



= W\j\- jdxj^x)A,(x), (38) 

A (x) - S -^& (39) 

Mx) ~ w (39) 

Taking the derivative of (38) over A it follows 
which after derivating over j v {y) and taking j M = A^ = leads to 



/ 



dz 



ST [A] 



S A^x)5A a (z) 



5 2 W[j] 



A=0 



-<W ^\x-y). 



(41) 



3=0 



Up to now the sources and field A^ have been considered as x^ depen- 
dent in eucliden space. Then, let us consider the static limit associated to 
equilibrium configurations. In that case 



3n 
A„ 



(42) 



(0,0,0,p(f)), 
(0,0,0, 0(f)), 

where it is assumed that when the external source is purely electrostatic, the 
system is able only to generate electrostatic mean fields. This is the case, 
in particular if parity invariance is not broken for example in the absence of 
external magnetic fields or other parity breaking effects. 
The effective potential is then given by the limit 



V[<£>(x)] = lim 



L 



r[o, o, o, 4> L {x , Xl )\ 



(43) 



in which <pL is taken as coinciding with <f)(x) in (0, L) and diminishing to zero 
outside. 

Now let us consider as in the previous section the variational problem of 
finding the extremum of the mean energy within the subspace of physical 
states under a restriction for a fixed mean field A^ = (0,0,0, 0(f)) . Due 
to the extremality condition a Schrodinger equation should follow for the 
ground state in presence of a source j M = (0, 0, 0, p{xj) in the form 
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(h + J d 3 x j,(x) \Qj) = EUWQj), (44) 

where as before H is the system Hamiltonian in the extended space. 
Taking the scalar product with \Qj > 

(n 3 \H\n 3 ) -E g = (E\j\ -E g - f d 3 x p(xMx)^j . (45) 

But in the euclidean limit of the Matsubara theory at great values of L 
it should follow 



- he in Z[j] ^ -he in 
W[0,0,0,p(x)\ = L E\j\. 



exp 



he 



(46) 



Then, in a similar way as in the previous section, and because \Qj > 
is a conditional extremal within the physical states subspace, the following 
inequality should be valid 



{VLj\H\VLj) - Eg = lim 



^ (V[0,0,0,p(£) - Jdx p{x)<p L {x, 



lim (\ r[0,0,0,y9 L (x,x 4 )]) -Eg 
V[(f{x)\ — Eg > 0. 



E n 



(47) 



Let us find now an expression for the kernel of the second derivatives of 
V which due to (47) should be a positive definite one. The Fourier transform 
of the equation (41 ) can be written as [8] 



' k 2 



(48) 



where IT^ is the polarization tensor which due to gauge invariance satisfies 
the tranversality condition 



U^k) = 

These relations allow to evaluate the Fourier tranform of the kernel 



(49) 



8 2 V\ 



8<t>(x)S<j>(y) ' 

Note that the increment of a functional under a X4 independent increment 
in the argument is given by 



8F[<t>\ 



00 5 pu] 

dx 4 - — - 5<f>(x). 

X, X4) 



(50) 



Thus, after using translation invariance 
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5 2 V[(j)\ 
5(t>(x)5(f){y) 



- [dx 4 dy 4 ^ 2r [°' ' '^ 



<f>=0 



dx 4 



5 2 rfo,o,o,^l 



<P=o 



8(j)(x — y, x 4 ) 50(0, 0) 

= r^(x-y,0), 
where r( 2 )(x — y, k 4 ) is the temporal Fourier transform of the kernel 

<$r[o,o,o,<? 



6<j>(x) 6<f>(y) 



(51) 



(52) 



4>=o 



But, integrating relation (41) over y 4 and taking the /j, — v — 4 component 
leads to 



J d 3 z rg(f-z,0) V aA (z-y,0) = -5^(x-y), 
where as before V^ v is the temporal Fourier transform of 

5 2 W 



$jn(x) 6j v (y) 
In momentum space it follows 



(53) 



(54) 



3=0 



r£ ) (fc,o)p a4 (fc,o) = -i. (55) 

However, from the Dyson equation (48) in the static limit k 4 = and 
/j, — v — 4 also follows 



fc 2 

— <5 4a - U 4a (k, 0) ) D a4 (A;, 0) = -1. 



(56) 



This relation can be further simplified by noticing that by transversality 
and the assumed parity invariance IL^ should have the structure [8] 



H-nv(k) = ix t (k) S,,,, - 



k^k u 



i/u)+TTu(k) V- 



(57) 



with 



Thus, in the static limit k 4 = 



(58) 



n 



M 




'-4H 



+ TX 44 {k)l 4 li 



(59) 
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Therefore, relation (56) reduces to 



47T 



-n 44 (/c,o) v u (k,o) = -1, 



which gives k 2 /An — n 44 (A;, 0) as the Fourier transform of the stability kernel 
T^(v — y,0) . Henceforth, the stability condition takes the form 



d 3 xd 3 y 4>{x) 
J d 3 k<f){k) 



5 2 V[(t)} 



5<f)(x)5(f)(y) 
k 2 



4ir 
k 



-n 44 (A;,0) <f>(k) 



J d 3 k 4>*{k) — (l - v(k)P(k)) (f)(k) 
P 

d 3 k <j>*{k) — e (k)(f)(k) > o, 

for arbitrary fields <f>(k) and where e(k) and P(k) are the dielectric and po- 
larization functions and v(k) = An/\k\ 2 . Relation (60) is the mentioned 
property. That is, the positivity of the dielectric function as an additional 
stability condition complementing the ones given in [3] and [4]. 



4 Summary 

It has been argued that the postivity of the dielectric function should be 
a necessary condition for the stability of a translational and parity invari- 
ant ground state of a charged many particle system at zero temperature. 
The connection of this condition with the alternative ones in [3] and [4] was 
analyzed. 

In complementing the here presented viewpoint it could be of interest to 
comment on the Wigner crystallization problem. It was the main example 
given in [3] and [4] for supporting the possibility for negative values of the 
dielectric function at zero temperature in a homogeneous state. The present 
discussion suggests that in lowering the charge density at T = the ho- 
mogeneity is broken at a point of infinite screening e— > oo which is the 
unset for the appearing of negative values of e(k). Then, if the positivity 
condition is correct, the translation invariance of the ground state is lost and 
the system pass to an inhomogeneous state (local order, plasma wave, ...). 
In such phases the notion of a dielectric constant with only one wavevector 

— * 

argument k lost the meaning. Then, in further diminishing the density, the 
system could attain the Wigner crystallization phase. Therefore, it seems 
that requiring as in [4], that after lowering the density, the Wigner crystal 
can only appear as related with a solution of a dispersion equation e(k) = 
for a homogeneous state ( only one k argument of the dielectric function) is 
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not necessarely valid. But, this condition is the main argument given in [4] 
in order to associate the Wigner crystallization problem as a physical system 
showing negative values of e(k) at zero temperature. 

Finally, the situation at finite temperature should be also commented. 
The present discussion and the previous work [5] indicate that a similar 
analysis could be valid at T 7^ 0. This problem will be considered elsewhere. 
In particular in Ref.[12], evaluations of the dielectric function at T 7^ for 
various classical plasma systems have been presented. The results show neg- 
ative values at inverse wavelength regions of the order of the interparticle 
distances. At such wavelengths the notion of a mean field as a statistically 
well defined thermodynamical internal parameter seems to loss sense. Then, 
under the assumption that a similar condition e(k) > could be found, as 
implied by a minimal value of a thermodynamical potential, the negative val- 
ues of the dielectric constant obtained in classical plasma systems could be 
associated with anomalous non-gaussian probabilities for fluctuations around 
the homogeneous field configuration [11]. It seems to me that such a situation 
could be interpreted perhaps, as a sort of weak breaking of the translational 
invariance possibly associated to some kind of measurable spatial inhomo- 
geneities. It should be noticed that the thermal energies related with the 
mentioned results are substantially lower than the coulomb energy at the 
mean interparticle distances. This observation leads to the interest in inves- 
tigating the local spatial order in molecular dynamic calculations for such 
systems ( the Monte Carlo configurations, not following the equation of mo- 
tion, may introduce an artificial homogeneization in the results). A further 
investigation of this more involved T^O case will be considered elesewhere. 
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